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Abstract 



A review of computations of free energy for Gibbs states on stationary but 
not static gravitational and gauge backgrounds is given. On these backgrounds 
wave equations for free fields are reduced to eigen-value problems which depend 
Q^' non-linearly on the spectral parameter. We present a method to deal with such 

^ • problems. In particular, we demonstrate how some results of the spectral theory of 

second order elliptic operators, such as heat kernel asymptotics, can be extended 
to a class of non-linear spectral problems. The method is used to trace down the 
^ ■ relation between the canonical definition of the free energy based on summation 

over the modes and the covariant definition given in Euclidean quantum gravity. As 
an application, high-temperature asymptotics of the free energy and of the thermal 
part of the stress-energy tensor in the presence of rotation are derived. We also 
discuss statistical mechanics in the presence of Killing horizons where canonical 
and Euclidean theories are related in a non-trivial way. 



These notes are based on the lectures delivered at the International Meeting on "Quantum 
Gravity and Spectral Geometry" in Naples, Italy (July 2- July 7, 2001). 



1 Introduction 

1.1 Basic facts 

Quantum field theory at finite temperatures appeared some forty years ago [|ll-@ and 
now it has many apphcations ranging from QCD to physics of the early universe, see 
|^-[0 for a review. These notes concern one particular aspect of this theory, namely, 
computation of the effective action and free energy in external gravitational and gauge 
fields. 

A finite-temperature field theory is closely related to Euclidean theory. This fact is 
now well known and is exploited in numerous applications ||^-[§[- To begin with we recall 
how to see this relation in a simple way. Consider a model of a real scalar field in 
Minkowski space-time with the action^ 



^t[0] = 1 t dt' f Sx [-(9t0)2 + (9,0)2 + m202 
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We study a thermal equilibrium of this field in a finite volume at non-zero temperature 
/3~^, a Gibbs state. The state is determined by the partition function 

Z{(3) = Tr e-^^, (1.2) 

where if is a normally ordered Hamiltonian of the given model. Normal ordering guar- 
antees that Z{[3) is well-defined. In general, ( p.. 2] ) also requires H to be an operator 
bounded from below (which is true for ( |1 . 1| ) ) . The average {(^\e~^^\(^) in some quantum 
state is equivalent to a transition amplitude from this state to itself for the period of time 
t = —if]. Such analytical continuation from the real to pure imaginary time is called the 
Wick rotation. One can define a complete set of states \ip) where the field operator is 
diagonal and takes the values ip. Then matrix elements of the evolution operator can be 
represented by a path integral [^ 

(yp|e-^*^|y.)=M|D0e-^^*['^], (1.3) 

where A/i is a normalization coefficient and trajectories (in a space of functions) in ( |1.3| ) 
begin and end at = y). As a result of the Wick rotation, It[4>] transforms into the 
functional 

/|[0] = z/_,40] = ^l^^'rfr I rf3^ [(9.0)2 + (9,0)2 + mV] • (1.4) 

The latter is the action for the same model but in a flat space-time where the metric is 
positive definite. One says that such as a space-time has the signature (+ + ++) and calls 
it Euclidean, as opposed to the physical space-time which has the signature ( — h ++) 



^We work in the system of units where h = c — 1. Our conventions for the curvature and metric 
coincide with the book [Q. 



and is called Lorentzian. Now to get ( p..2| ) from ( |1.3| ) one has to take the trace over all 
possible values 93. Thus, the partition function can be written as 

Tx e-^^ = mJ D(j) e-'f^^K (1.5) 



By construction, the path integral is taken over all possible closed trajectories going along 
the Euclidean time r. This is equivalent to imposing periodicity on r with the period 
13. This result can be easily extended to non-zero spin fields including fields with Fermi 
statistics. 

The integral in r.h.s. of ( |1.5|) is exponentially suppressed for large values of and this 
allows its rigorous definition. One can also introduce interaction of the field with an 
external source. Then the Euclidean path integral can be used as a generating functional 
for Euclidean Green functions 1^,0. The Euclidean Green functions correspond to real- 
time finite-temperature Green functions which obey certain boundary conditions known 
as Kubo-Martin-Schwinger (KMS) condition 0, see [|l-[§[ for more details. 



It is instructive to give another derivation of (1^). Let us introduce the free energy 
of the system 

F(/?) = -/?-MnZ(/3). (1.6) 

To calculate it consider a single-particle excitation (puiit, x) = e~*'^*0aj(x) of the field with 
a frequency u. The spectrum of u can be found when (puj(t, x) is substituted in the wave- 
equation (in Lorentzian space-time). As is easy to see, it yields an eigen- value problem 

^'0. = uj^<Pu, (1.7) 

H^ = -d'f + m\ (1.8) 

The operator H = VH"^ is a relativistic analog of the quantum-mechanical Hamiltonian. 
In the considered case of a "system in a box" the spectrum is discrete and positive. F{i3) 
is well-defined and can be represented in the form 

F{f3)=p-'Y.^n{l-e-^^), (1.9) 

which is true for free fields. 

Let us also define the Euclidean effective action 

W^{(3) = -In f D(l) e-^/^['^l = -^InA , (1.10) 

J 2 ^ 

where A are eigen- values of the Euclidean operator L^ = ~d^ — df + rn?. The series in 
( |1.10| ) diverges and as usually has to be regularized according to some prescription. The 



eigen- vectors of Le are periodic in r. They have the form e^^^'^cf)^ where o"/ = (27r/)//5 and 
/ is an integer. Hence, A = af + w"^ where uo are defined by (|1.7] ). The quantities ai are 
known in the literature as the Matsubara frequencies [ffl. 



The relation between F{f3) and W (jS) can be found by using the identity |]TO[,|n 



In (l - e~^'^) = 



- hm —-Ciulu, p) -— 



2 ' uj' 
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(1.12) 



l=—oo 



The C-function ( 1.12 ) is defined at z/ = by analytical continuation from the region 



Re u > 1/2. By taking into account ( |1.9|) we find that 



pF(p) = lim 

UIQ— >00 



E 



W<W(i _ 



— hm -—C{i'\uj, p) 



;i.i3) 



The first term in r.h.s. of this equation formally coincides with the determinant of the 
Euclidean operator, see ( p..lOD . The second term in ( |1.13| ) is related to the energy of 
vacuum fiuctuations 



Eq 



E 



u 



Hence we conclude that 



(3F{(3) = W'^{(3)-(3Eo 



(1.14) 



:i.l5) 



Although both W^ and Eq diverge their difference is finite and coincides with the free 
energy. In what follows we give more rigorous definition of W^ and Eq. Obviously, ( 1.15 ) 
is in agreement with relation ( |1.5| ) derived by the path- integral method where In A//? is 
associated to [3Eq. 

In case of constant background fields the Euclidean effective action is closely con- 
nected to the notion of effective potential |[T^ and is a very useful tool in studying phase 
transitions in a system. Relation ( |1.15|) can be extended to situations when the system is 
placed in an arbitrary static gravitational field as was first done in pioneering works by 
Gibbons |jl3| and Dowker and Kennedy [|1^ (see also [jl5[ for discussion of ( p..l5| )). This 
opens the possibility to study a variety of physical phenomena in gravitational fields. In 
a flat space-time the vacuum energy Eq depends on the boundary conditions, typically 
on the size of the system. In more general situations £^0; like other thermodynamical 
characteristics, becomes a functional of external background fields. Thus, the difference 
of W^ and F{[3) is essential. 

We discussed a finite-size system. In many physical situations, however, one has 
to consider the thermodynamical limit when the size becomes infinitely large. Strictly 
speaking, in this situation the Gibbs state cannot be realized by a density matrix in a Fock 
representation. The partition function is not defined because the operator exp(— /?if) is 



not of trace class. As was shown by Haag, Hugenholtz and Winnik |]16| the mathematically 
satisfactory definition of a Gibbs state in this case can be based on the algebraic Gel'fand- 
Naimark-Segal (GNS) construction. The state is realized as a reducible vector state in 
GNS representation compatible with KMS condition. 



Although the density matrix cannot be introduced for spatially infinite systems the 
free energy (|1.9|) and some other quantities like Green functions can still be defined. One 
can start with a "system in a box" and consider a suitable limit when the box is taken 
away to infinity. Because of spatial homogeneity free energy of the considered model will 
be proportional to the volume of the box. Thus, its density per unit volume is well-defined. 
Analogously one can define the density of F[/3) when the spatial part of space-time is a 
hyperbolic manifold of constant negative curvature. This can be done by replacing the 
summation over the modes in (|1.9| ) with an appropriate integral where the integration 
measure is determined by the spectral density of the operator H^, see [|17| for a review. 
We will follow a similar method in next sections. 



Equation ( |1.15| ) will be a starting point of our further discussion. It relates the two 
definitions of free energy: the first is canonical, F^ {(3) = F{(3). It is based on (|L^) and 
is called "summation over the modes". It is the definition which follows from principles 
of statistical mechanics and has a clear physical meaning. The second is the definition of 
free energy in terms of the Euchdean effective action (|1.1CI| ), i.e., as F^{(3) = P~^W^{P). 



The two definitions are related^] as (see ( |1.15| )) 

F^(/5)=F^(/5) + Eo. (1.16) 

Once restricted to static space-times W^{P) and, hence, F^[[3) have a number of attrac- 
tive properties, i) The classical action I§ has a covariant form in Euclidean space-time. 
Under a proper choice of integration measure in the functional integral the effective action 
becomes a covariant functional of the background metric, the property which is especially 
important in quantum gravity, ii) Euclidean wave operators L^ are positive elliptic op- 
erators. There is a number of rigorous mathematical results concerning the properties 
of these operators |T3] which can be used to define F^ . iii) Contribution of the vacuum 



energy Eq is included in F^ . This is one of the reasons why F^ , as distinct from F*-^, can 
be defined as a covariant functional. In principle, like in a classical theory, by varying F^ 
over the background metric one can define the total finite-temperature stress-energy of a 
quantum field. 

1.2 Difficulties 

We now discuss situations where the Euclidean definition of free-energy faces difficulties. 
Consider the same model as before and assume that the field is inside a cylinder of radius 
R. If the cylinder has a finite length we still deal with the field in a box. Let us define the 
Gibbs state as a thermal equilibrium of the field as measured in the frame of reference 
which rigidly rotates around the symmetry axes of the cylinder with the angular velocity 
VL < R~^. The metric in the rotating frame can be obtained from the Minkowski metric 



^In some works F^ is also defined witliout the vacuum energy Eq. In this case on static space-times 
F^ and F^ coincide. 



written in polar coordinates (r, ip) (with the origin at the center of the cyhnder) where 
the polar angle ip should be replaced to (p = ip — fit 

ds^ = -(1 - r^Q^)dt'^ + 2Qr^dtd<p + r^d<p^ + dr'^ + dz^ . (1.17) 

The free energy can be represented in the form (|1.9|) where summation is taken over 



frequencies ui of quanta with the angular momentum along the z axis equal to / (/ is an 
integer). If uji is a frequency in the non- rotating frame, then uJi = uJi + IQ. This relation is 
easy to see by rewriting a single-particle solution with momentum / in coordinates ( |1.17| ) 



The partition function of the system is defined by (|1.2| ) where H is & normally ordered 
Hamiltonian operator generating time evolution in the rotating frame ( |1.17|) . It is not 
difficult to construct a formal Feynman path integral representation for this partition 
function and come to equation (pTSJ ). To get the Euclidean action I^ one has to start 
from action It in Minkowsky space-time, write It in coordinates ( |1.171 ), go to Euclidean 
time T = it and impose periodicity on r with period j3. One finds that 



/|[0] = I^,p[(t)] = J dr J r^drd^dzlg^-'d^cPd.cP + m^'] , (1-19) 



r/3 
/o 
where g'^'^ is a contravariant tensor defined as 



g^'p^.Pu =pI + 2tnprp^ + (r'2 _ ^2)^2 +pI+pI_ (1.20) 

The important difference as compared to the non-rotating ensemble is that the form 
( |1.2QD is complex. Thus, the Euclidean classical action in the Feynman path integral (|1.2CI| ) 



becomes a complex functional. This results in a difficulty because the corresponding wave 
operator in the path integral is not a positive-definite elliptic operator]^. The Euclidean 
effective action W^ for such complex operators cannot be rigorously defined and the 
Feynman path integral (|1.5|) is nothing else but a formal expression. 

This problem can be avoided if we simply make an additional step and consider analyt- 
ical continuation to pure imaginary values oiQ = iCl. Then the Euclidean action becomes 
real, the Feynman integral converges, and the Euclidean operator Le is a positive-definite 
elliptic operator with the leading symbol 

r''P^.P. =pI- ^np^P^ + (r-2 + fi2)p2 + p2 ^ p2 (^21) 

Now we can interpret our theory as truly Euclidean, i.e. on a real space-time with the 
signature (+ + ++) and metric 

dsl = g^^dx^dx"" = (1 + r'^f?)dT'^ - 2Qr^dTd<p + r^d<p^ + dr^ + dz^ . (1.22) 



^Positive-definite elliptic operator is an operator whose leading symbol ( 1.2C ) is a positive-definite 
quadratic form for all p ^ 0, see [nSl. 



By going to the new coordinate ip = (p — VLt one can check that ( |1.22| ) is flat. Thus, the 
Euchdean manifold is a flat R'^ space where points (r, ip, r, z) and (r + /5, yj — il/3, r, z) are 
identified. 

This simple flat space example is a good illustration of what happens when one is trying 
to relate statistical mechanics on an arbitrary stationary space-time to the Euclidean 
theory. Stationary but not static metrics look as 

ds^ = gttdt^ + 2gtidtdxi + gijdx^dx^ (1-23) 

where components g^y do not depend on t. Such metrics appear in different physical prob- 
lems. In particular some black hole or cosmological solutions of the Einstein equations 
(e.g., Kerr, Kerr-Newman or Godel metrics |T^) have form (|1.23|) . The Wick rotation in 



this case should be accompanied by an additional analytical continuation of the param- 
eters of the metric to get from (|1.23|) a real Euclidean metric. This procedure implies 
complexification of the space-time, the idea which was first suggested and used by Hartle 
and Hawking ||2^ and Gibbons and Hawking [^|,|22]. More precisely, this procedure 



can be described as follows [^. One considers a 4- complex-dimensional manifold A^c 
with a complex contravariant tensor field g'^^ of type (2,0). Lorentzian and Euclidean 
manifolds are defined as real sections of A^c, i-e., as 2-complex-dimensional submanifolds 
(slices) on which the restriction of g'^'^ to the cotangent space of the slice is a real form 
with the signature ( — h ++) or (+ + ++), respectively. It should be noted that an ar- 
bitrary Lorentzian space-time may not admit a complexification which has a real section 
with Euclidean signature and vice versa. We will exclude these cases from our further 
discussion. 

Although, the described procedure of going to Euclidean theory makes the Feynman 
integral well-behaved one loses the original relation to statistical mechanics. The problem, 
thus, is to find out whether there is an analog of equation ( |1.15|) when the space-time is 
stationary but not static. 

There is another problem which is not present in static space-times. In static space- 
times the spectrum of frequencies u is determined by equation ( |1.7|) . The latter is obtained 
from the wave equation for a field under the substitution 0w(t, x*) = e~^'^^(f)i^{x). Equation 
( |1.7|) is the eigen-value problem for H^ which is a 3-dimensional elliptic operator (see 
( |1.8|) ). However, on a stationary space-time the same substitution results in equation 

(-^"(x)cu2 + iuj{-2g''{x)di + b{x)) + g'^{x)d,dj + c(x)) 0^(x) = 0, (1.24) 

where b{x) and c(x) are defined by the metric and parameters of the model. This is an 
eigen-value problem which depends polinomially on the spectral parameter u. Equation 
( |1.24| ) belongs to a class of non-linear spectral problems which appear in different physical 



situations. A typical example of such a problem is a two-particle scattering where the 



scattering potential depends on the energy of particles [0]. The difficulty with ( |1.24| ) is 



that its quantum-mechanical interpretation is obscure but, what is worse, this equation 
is not an eigen-value problem of any operator. Thus, standard results from the theory of 
elliptic operators cannot be used here. Certainly, it is possible to analyze ( |1.24| ) in each 
particular case and try to solve it explicitly. This is not easy, however^. For instance, in 
case of rotating black holes ( p..24| ) is rather complicated and the fact it allows separation 
of variables p^ is a big luck. For fields with spins 1/2, 1 and 2 the situation gets more 



complicated, see [p^] . 

Similar difficulties occur in quantum theory in external static gauge fields. In the 
presence of electric field the time component At{x) of vector potential is not zero and 
time derivatives in wave equations are covariant derivatives dt + ieAf. To have after the 
Wick rotation a positive definite elliptic operator one has to go to imaginary potential 
At = lAr- This requires complexification of the vector field A^. Analogously, like (|1.24|) , 
the equation for spectrum of frequencies u of single-particle states in this situation is not 
a standard eigen-value problem because it contains terms linear in uj. 

Finally, we dwell on another difficulty which appears in case of black holes. A remark- 



able property of black holes is the Hawking effect p9|: in quantum theory black holes 
evaporate by emitting particles with a thermal spectrum. The corresponding temperature 
(the Hawking temperature) is T^ = k/{2it) where k is a characteristic of the gravitational 
field strength near the black hole horizon called the surface gravity (for an introduction 
to black hole physics see 0],|2O|)- '^^^ Hawking effect implies that a black hole can be in 



thermal equilibrium with a quantum gas if the temperature of the gas coincides with the 
Hawking temperature (for a Schwarzschild black hole the corresponding state is known 
as a Hartle-Hawking vacuum pO|). 



This result can be understood on purely geometrical ground. The Schwarzschild metric 
is a static spherically symmetric metric 

ds^ = -B{r)dt^ + — -- + r^(sin2 Od^'^ + dd'^), (1.25) 

B[r) 

where B{r) = 1 — Vh/r. The horizon is located at r = r/^ and is a 2-sphere of radius 
r/i = 2MG, where M is the mass of a black hole, G is the Newton constant. The 
Euclidean theory corresponding to a finite-temperature theory around a Schwarzschild 
black hole is defined on a space obtained from ( |1.25|) by changing t to —ir where r has 



to be periodic with a period (3. As is easy to see, near the horizon the Euclidean metric 
looks as 

dsl = K^p^dr'^ + dp^ + rlisi-n^ Od^/ + dO"^), (1.26) 



where p ~ 2Jrh{r — r^) and k = l/(2r/i) is the surface gravity. The coordinates p, kt in 
( |1.26|) behave as polar coordinates on the plane. If f3 is arbitrary the space with metric 
There is a number of publications in mathematical literature where spectrum and eigen functions are 
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studied for one-dimensional non-linear spectral problems similar in the form to (1.24), see, e.g., p5[ ,| 
and references in fE6|. 



K^p^dr'^ + dp^ is a cone. The conical singularity vanishes at /5 = 271/ n. This corresponds 
exactly to the Hawking temperature. 

\i [3 ^ Ph = T^^ the Euclidean space is not smooth. Conical singularities are defects 
of the geometry where the curvature has delta-function-like behaviour. In the presence of 



such defects quantum theory acquires new ultraviolet divergencies |£T|-|3^ in a form of 

additional divergent terms in the effective action with a support on the Euclidean horizon. 

The canonical formulation of statistical-mechanics in the presence of Killing horizons 

also has non-trivial features. The local temperature measured at a point with coordinate 



r is given by the Tolman formula, T{r) = Th/ J B{r). The factor yB appears because of 
the gravitational blue shift of frequencies near the horizon with respect to frequencies at 
infinity. When r approaches r^ the local temperature becomes infinitely large. So near 
a black hole one has a thermal bath at very large temperature. In this regime masses of 
fields play no role and the canonical free energy F'~^ acquires infrared divergencies. 

As a result, in the presence of horizons the Euclidean and the canonical free energies 
look different and finding the relation between them becomes a problem. 

1.3 Content 

The remaining part of these notes is mainly devoted to two topics where the results are 
new, i.e. to studying non-linear spectral problems and finding relation between canonical 
and Euclidean methods for stationary backgrounds. In Section 2, we present an approach 
which enables one to reduce non-linear spectral problems to standard eigen- value problems 
of some fiducial elliptic operators. In particular, in section 2.3 we demonstrate that some 
elements of the spectral geometry can be extended to this class of spectral problems. We 
introduce an analog of the trace of the heat kernel operator (a pseudo-trace) and show that 
its asymptotic expansion is determined by local geometrical invariants which generalize 
the known heat-kernel coefficients. We briefly discuss some physical applications such as 
high-temperature asymptotics of the free energy and of the stress-energy tensor including 
the case of non-zero rotation. In Section 3, we use our method to study relation between 
statistical mechanics and Euclidean theory on arbitrary stationary space-times. Among 
other advantages, the Euclidean theory allows for an alternative interpretation of the 
results of Section 2. In Section 4 we discuss possibility to extend the spectral asymptotics 
to a larger class of non-linear spectral problems. In particular, we demonstrate how it 
can be done in case of an external electric fleld. The aim of Section 5 is to consider 
the features which appear in case of quantum theory in the black hole exterior, or, more 
exactly, in case of Killing horizons. 

It should be emphasized that we are not trying to present mathematically rigorous 
proofs of our statements but rather use arguments which are closer to physicists. Thus, 
some of our conclusions can be considered as conjectures. We test, however, how they 
work in concrete non-trivial examples. Our discussion will be restricted to case of free 



scalar fields. This enables us to avoid unnecessary complications and concentrate on the 
main problems. When possible we give references to results concerning higher spin fields. 

2 Stationary space-times 

2.1 Killing frame 

We begin with some helpful definitions. Consider a field on a domain Ai oi a. D- 
dimensional space-time with a time-like Killing vector^ field ^^ (^^ < 0). Ai may be a 
complete manifold if ^'^ is everywhere time-like. In most other cases C,^ is time-like only 
in some region. We will study solutions of field equations in the frame of reference related 
to Killing observers whose velocity m^ is parallel to ^^ 

u^ = B-^'^e , B = -^\ (2.1) 

For a Killing observer a solution 0^^ carrying the energy u is defined as 

iC^cpu; = ujcp^ (2.2) 

where C^ is the Lie derivative along ^^. The background metric (7^1, can be represented as 

gf,u = h^^ - u^u^, (2.3) 

where h^^ is the projector on the directions orthogonal to u^. Because the shear and 
expansion of the family of Killing trajectories vanish identically the trajectories are char- 
acterized at each point only by their acceleration w^ and the rotation A^i, with respect 
to a local Lorentz frame 

w^ = u^.xu^ , (2.4) 

V = 2^i^'("^;p-"p;^)- (2.5) 

To proceed it is convenient to choose coordinates x^ = {t, x*), i = 1,D — 1, where ^ = d/dt 
and, consequently, ho^ = 0. According with (|2.1|) , ( pTSf ), the interval on Ai can be written 
as 

ds"^ = -B{dt + ttidx'f + dt^ , (2.6) 

'Ufj,dx^) = dt + ttidx^ , (2.7) 

(2.8) 



B 

u 



'B 

dt^ = h^^dx'^dx" = hijdx'dx^ . (2.9) 



'By definition the Killing vector is a solution of equation S,p,-i, + ^i/;^ = 0. 
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The role of Oj becomes clear under synchronization of clocks by sending light signals: the 
two events with coordinates (t, x*) and (t — akdx'',x'' + dx*) are occuring at the same 
moment. The metric dP serves to measure the proper distance between these points. 

In coordinates (t, x*) the only non-zero components of acceleration ( p.4|) and rotation 
(|3D are 

Wi = -{\nB)^i , Aij = ~-y/B{aij - a^^i). (2.10) 

In four-dimensional space-time one can define a vector of local angular velocity 

a = ]^e,,kA'\ (2.11) 

where eijk is a totally antisymmetric tensor. The absolute value of the angular velocity is 

fi = (afi*)i/2 = Qa^-A^,) . (2.12) 



The form of the metric in the Killing frame, equations ( |2.6| ), (|2.7] ), is preserved under 
arbitrary change of coordinates x* provided hij and Oj transform as a D — 1 dimensional 
tensor and a vector. There is also another group of transformations, which preserves 
( p.6|) , ( p.7|) , namely, t = f + f{x), ai = a[ — dif{x), where / is an arbitrary function of x\ 
Under these transformations a, changes as an Abelian gauge vector field. By considering 
single-particle excitations with fixed energy u 

(j)Ut,x') = e~"^'Mx') , (2.13) 

one can realize this group of transformations as a local U{1), 

<jy.{t,x')=<Pi{t',x')=e-'-UUx') , 

U^') = e'-f^^UUx') . (2.14) 

In this picture, uj coincides with an "elementary charge". To quantize in the Killing frame 
one needs a full set of modes 0i^(x*). As follows from the above arguments, the equations 
which determine (pujix^) have a form oi D — 1 dimensional equations for charged fields in 
external gauge field Oj on a space with the metric hij. It is important that covariant prop- 
erties of the theory in D dimensions guarantee diffeo- and gauge-covariant form of the D — 1 
dimensional problem. Such a reduction from Z) to Z) — 1 is analogous to the Kaluza-Klein 
procedure which yields the Einstein-Maxwell theory from higher dimensional gravity. The 
difference between the two reductions is that in the standard Kaluza-Klein approach the 
"extra" dimensions are compact and the charges are quantized. 

Let i3 be a Z^ — 1 dimensional space with metric dP, see ( p.9|) . Consider a point p on 
B with coordinates x* and a vector Vi from the tangent space at p. On Ai, p corresponds 
to a trajectory of a Killing observer with the same coordinates x*. At any point of the 
trajectory one can define a vector V^ orthogonal to u^ such as V^ = /if V^. Suppose that 
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connection V^ on B is determined by hij. Then the covariant derivative with respect to 
this connection can be written as [^ 

V,V, = hfh''^Vx,, , (2.15) 

where V^.^, is the covariant derivative on A^ with respect to the connection defined by g^j,^. 
(One can easily check that Vkhij = 0.) Relation ( ^.151 ) can be generalized to an arbitrary 
field on Ai. For instance, for a scalar field 

h'^d^cj) = {dj - a^dt)(j) = D4 , (2.16) 

for a vector orthogonal to u^ 

h'>h';V^,p = (V,- - ajdt)V., = DjV, , (2.17) 



where Vi = hfV^. The time derivative in ( p.l6| ), ( p.l7|) appears in general because fields on 



Ai change along the Killing trajectory. If and V^ are solutions with certain frequency, 
see ( p.l3|) , then Di become covariant derivatives on B in external gauge field Oj. This 



demonstrates explicitly diffeo- and gauge-covariance of the theory which are left after the 
reduction. 

2.2 Kaluza-Klein method 



We now present an approach developed by Fursaev in |3^ Q. Consider a real scalar field 
(j) which satisfies the equation 

i-V^V^ + V)(j) = 0, (2.18) 

where \^ is a potential. Suppose that A^ is a globally hyperbolic space-tim^ and S is a 
Cauchy hypersurface in Ai. One can define the Klein-Gordon inner product on E (see 



< 01,02 >=/^f/S%(0i, 02), (2.19) 

Jm(</'i. <t>2) = -^{<Pl^^<p2 - 9^0t02). (2.20) 



If 01 and 02 are solutions to ( p^.lSD the current j^ is divergence free, V^j^ = 0. This 
property guarantees that ( |2.19D is independent of the choice of E. 



Consider first systems in a finite volume (space B is compact). Then the problem 
is formulated as follows: one has to find a complete set of solutions (f)^(t,x^) to ( |2.18| ) 
with certain frequencies u which are ortho-normalized with respect to the product ( [2.191 ). 

^ We discuss scalar fields only. Ref. Q also includes the analysis of the spinor fields but has a mistake 
in the definition of a spinor Hamiltonian. Thus its results for spinor fields require a revision. 

^A globally hyperbolic space-time is a one which admits a Cauchy surface. A Cauchy surface E is a 
space-like hypersurface such that every non-space-like curve intersects S exactly once. For more details 
see |l§. 
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Frequencies uj are energies of single-particle excitations. The spectrum of uj can be used 
to define the canonical free energy F^ in the Gibbs state by ( |1.9|) . 

Let us demonstrate how finding the spectrum can be reduced to a standard eigen- value 
problem for an elliptic operator. By using the results of the preceding section the wave 
operator in the Killing frame (p.3|) can be represented as 

V^V^ = g^'D^D,, (2.21) 

^^ = V^ - a^du (2.22) 

where a^dx^ = aidx\ The connections V^ are determined on some space Ai with metric 

ds'^ = g^^dx'^dx'' = ~Bdt^ + df. (2.23) 

Relation between M. and M. becomes transparent when comparing ( p.23|) with ( p.6| ) . We 
will call M. and a^ the fiducial space-time and the fiducial gauge potential, respectively. 
Let us consider now a scalar field ^'■'^^ on tVI which obeys the equation 

(-^'"(V^ + i\a^)iy, + i\a,) + V)0(^) = 0, (2.24) 



where A is a real parameter. Because ( |2.24| ) does not contain terms linear in the time 
derivatives, it takes the following form 

H\\)cl^^^\x^)=uj\\W{x^) (2.25) 

on functions 0^^^(t,x*) = e~*'^*0£^''(x*). The operator if (A) has the meaning of a relativistic 
Hamiltonian for the field (j)\^\x^) on B. It takes the simplest form after the transformation 
(which does not change the spectrum) 

if (A) = e-^"ii(A)e^", (2.26) 

H'^[\) = -WiVi + iXai){yj + iXaj) + V, (2.27) 

e-"^" = -^2 = B. (2.28) 

Connections Vj correspond to fields on a L* — 1 dimensional space B conformally related 
toB 

df = hijdx'dx^ = e^^df. (2.29) 

The "potential term" in (2.27) is 



V = B 



D — 2 D — 2 



(2.30) 



where w^ is the acceleration ( |2.4|) . One can arrive at ( 2.27 ) by doing a conformal trans 



formation in equations (p.l8|) , ( p.24| ). Under this transformation the physical metric gfj_u 
changes to g^^, = g^u/B and the Killing vector on the rescaled space has the unit norm, 

e = -1- 
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By taking into account ( |2.27| ) we conclude that ( |2.25| ) is the standard eigen-value 
problem in a finite volume for a positive-definite elliptic operator. The spectrum of H^{X) 



is discrete and bounded from below |jT8[. For simplicity we assume that the spectrum is 
strictly positive. 

Because the fiducial spectrum depends on A we can find the physical spectrum from 
the constraint 

x{X,uj) = 0, A>0, (2.31) 

xiX,uj) = -uj\X)+X\ (2.32) 

If A obeys ( p.31|) we denote the corresponding solution to (|2.24| ) as (p^\t, x*) = e~*'^*0^^(x*) 
As follows from (p.21|) - (p.25|) , (p.31|) this field is also a solution to physical equation (|2.18|) , 
i.e. one can write 

(l)^{t,x')=C^(l)^fit,x'), (2.33) 

where C^ is a normalization coefficient. Therefore, the relativistic eigen-value problem 
can be reduced to the standard eigen-value problem for a one-parameter family of elliptic 
operators H^{X) and to solving (|2.31|) . We call this method of finding solutions to (|2.18| ) 
the Kaluza-Klein (KK) method. 

To fix C^ we have to analyze the inner products. The fiducial fields obey equations 



( p.24| ) which dictate a different form of the corresponding vector current and the product 



jf,{(pi, 02) = -i{(t)*i{d^ + iXa^)(t)2 - (<9^ - iXa^)(j)l(t) 



fM)ViHJ2), 



n,(P2) 



dt^~3,{ 



n,(P2) 



(2.34) 
(2.35) 



The connection between ( |2.19| ) and ( |2.35| ) can be established as follows. (Without loss of 
generality we put B = — ^^ = 1 because one can always reduce the problem to this form 
by conformal rescaling.) Choose S and S as constant time hyper-surfaces0. Then 



< 0c^,0a > = 



Vhd 



D-l 



X 



{u + a)(j)*^(j)a + i(f)la\Vi + iaai)(j)„ - i{Vi - iuai 



h*a\ 



W,0W) = (^ + a)/v^d^-y0W)>W, 



(2.36) 
(2.37) 



where h = det hij, and hij is the metric on B. Indices i,j are raised with the help of h'^K 
Note that up to a multiplier ( ^.37| ) coincides with standard product on space of square 
integrable functions L^. Equations ( p.36| ), ( p.37| ) can be used to write 



< 4>iu, 4>a >= {4>u,, 4>a) 

When ijj ^ o one can use the identity 

E'^{y)) - H^{a) = (to - a)(-2m*Vi - i^ia' + (to + a)a'ai 
®For fiducial space-time Ai metric induced on E coincides with metric on B. 



(2.38) 



(2.39) 
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to rewrite ( p.38| ) as 

1 

+ 



< 



^a>-- 



2(a^-uj^) 



(.., (H^co) - H\a))<P^y + (0^, {H\uj) - H\a)W. 



(2.40) 



The r.h.s. of ( p.40| ) vanishes if -ff^(A) is Hermitean, which we assume to be the case. ( 2.401 ) 
demonstrates that physical modes 0^^, (p^ obtained from fiducial modes under restriction 
( |2.31|) are orthogonal with respect to ( p.36|) when uo ^ a. 

Consider now the modes with equal frequencies, (p^ and 0^^. ( p.38[) can be written as 






By differentiating both sides of ( p.25| ) with respect to A one finds 

a,if2(A)0W + //^(A)9,0W = 9,^^(A)0W + ^^(A)9,0W, 

dxH^{X) = -2ia*(V, + iXtti) - iV'ai. 
Relation ( p.41|) can be rewritten by using ( p.42| ), ( p.43| ) and definition (|2.31|) as 



< 0a;, 0c^ > = 



2uj 



■^LO") y^uj J -J 



(2.41) 

(2.42) 
(2.43) 

(2.44) 



where x'(ti;,u;) = 9ax(A,c<j) at x(A,ct;) = 0. Because the form {4>^,(p^) is positive-definite 
it follows from ( |2.44] ) that the form determined by the Klein-Gordon product is positive- 
definite if 

X'iuj,iu)>0. (2.45) 

If x'(a;,u;) < for some u the corresponding state has negative Klein-Gordon norm and 
should be excluded. We will assume that condition ( p.45| ) is satisfied for any solution of 
( p.31| ) . We will see later that ( |2.45| ) also plays a role in the Euclidean formulation of the 
theory. 

If ( |2.45|) is satisfied then the orthogonality of (p^, 0^ with respect to the standard prod- 
uct implies their orthogonality with respect to ( p.44| ). Finally, if 0^-* have a unit norm, 
then according to ( |2.44| ), normalization of 0^^ in ( |2.33| ) requires |C^|~^ = x' {^ ^ ^) / C^^) ■ 

To summarize: we have proved that functions (|2.33| ) are a set of modes orthonor- 
malized with respect to the Klein-Gordon inner product provided the fiducial modes are 
orthonormalized with respect to ( p.35|) . This fact also means that there is no contra- 
diction between orthonormalization procedures for physical and fiducial modes. Given 
this result, one can go further and study other propertiers of the set of physical modes 
including its completeness 0. 

^For one-dimensional non-linear spectral problems (for Sturm-Liouville operator) the completeness of 
eigen- functions is discussed in |2q|. 
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2.3 Spectral coefficients 

We discuss now properties of the spectrum of single-particle states. For the class of models 
considered in section 2.2 one can define the function 

Kit) = Y.e~'-\ (2.46) 



where t > 0. If the space-time is static ( p.46|) coincides with the trace of the heat kernel 
of the operator H^, 

K{t) = TT e-^"\ (2.47) 



In general, however, K{t) is not a trace of any operator like (|2.47|) . Let us call it a pseudo 



trace. For an elliptic operator H^ the asymptotic of (|2.47 ) at small t has the known form. 



To avoid unnecessary complications we assume that the system is defined on a compact 
space without boundaries. Then at small t 

1 oo 

^ ^'"" ~ (4.t)(°-W2 g/"'"- (2,48) 

The coefficients a„ are called Hadamard-Minackshisundaram-DeWitt-Seeley coefficients, 
or sometimes the heat coefficients. They are local invariant functionals polynomial in 
curvatures of the background metric. 

The expansion ( p. 481 ) plays an important role in mathematical [Q and physical appli- 



cations. It is natural to ask whether there is an analog of this expansion for pseudo-trace 
K{t) in general case. 

To answer this question let us note that K{t) diverges at small t as a result of sum- 
mation over arbitrary large u"^. When t is finite there is the exponential cutoff of the 
series at cj^ ~ 1/t. Therefore, asymptotic properties of K{t) are related to the behaviour 
of the spectrum at large uj"^, the same property which determines the asymptotics of the 
heat kernel in ( p.48| ). For large uj the spectrum becomes sufficiently dense and can be 



considered continuous. We will use this fact for our purpose. Define K{t) for operators 
with continuous spectrum as 



K{t) = / <l>(cu)rfcue-*^ . (2.49) 

Here /z is the mass gap of the spectrum. To avoid complications we also assume that /i > 
and there are no isolated discrete levels[^ The quantity $(co') is the spectral density or 
the total spectral measure defined on a compact domain E^ of a Cauchy hypersurface S 
as 

$(w) = J^ rfS^'^j^(0^,;,0^,;), (2.50) 
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These assumptions seem to be not essential in studying the limit of large u>. 
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where j^ is defined in ( |2.20| ). We assume that (p^^i is a complete set of single-particle 



solutions to ( 2.18 ). Index I enumerates all modes with the same frequencies. In problems 
with a continuous spectrum E has an infinite volume and normalization of modes is to be 
understood in terms of distributions, i.e., as 



< 



^uj,ki 



^a,i >= 5ik5(u -a). 



(2.51) 



where 5{x) is the delta-function. The integral on r.h.s. of (|2.5CI| ) is divergent. This is the 
divergence of the large volume or infrared divergence. To avoid the divergence one has to 
work with a regularized density obtained by restricting the integration in ( p.50| ) to some 
compact domain S^ in S and then expanding S^ to S |^. 

The important advantage of the KK-method is that the spectral density ^{uo) can be 
related to the spectral density of the corresponding fiducial problem ( |2.25| ). Suppose that 



c(A) 



are a complete set of ortho-normalized solutions of (|2.25|) 



kW aW 



:i,<l>];j)=6i,6iuj-a). 



(2.52) 



If the modes are related by 

0^(t,xO=0(rHt,xO, (2.53) 

one can show Q] by using (|2.36| ) and ( |2.37|) that the normalization (|2.51| ) follows from 

(Pi. 

The (regularized) spectral density of fiducial modes is 



^u;;X)=J^ dJ:^Y.J^^i 



a<&), 



(2.54) 



where j^ is defined in ( p.34| ) and S^ is a compact domain in S. The definition of S^ 
corresponds to the definition of the domain S^ in S (we choose again S and S as constant 
time hypersurfaces) . Let us introduce also an auxiliary quantity 



vl/(^; A) = $(c.; A) - ^ E H'l d,H\X)<j>l'l) + {<P% d,H\X)<pl'l) 



(2.55) 



where dxH^{\) is given in ( |2.43| ). As follows from ( |2.38D 

$(cj) = ^(cjiu;). 



(2.56) 



It should be noted that $(c<j) does not coincide with ^{u; u), as one could naively expect. 
The distinction of these two quantities is in different forms of the inner products. Consider 
now the spectral representation for the heat kernel of H'^{X) 



Tr e-*^'(") 



^{uj; X)e~''^ duj 



(2.57) 



M 



-'^^We will see that this procedure enables one to study local characteristics of the system such as the 
density of the free energy per unit volume or a stress-energy tensor. 
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where integration in the trace is restricted by S^. We assume that the mass gap of H^{X) 
is positive and does not depend on A (which is true in a number of physical problems). 
Similarly we can define the integral 

1 



/ ^{uj- A)e-*'^ duj = Tr 



1 - ^d,H\X) ] e"*^^W 



(2.58) 



where the r.h.s. is the consequence of ( p.55|) . Because the trace does not depend on the 



choice of the basis and, hence, on A one can write ( |2.58|) as 



r ^{uj; X)e~'^'duj = (1 + ^dx) Tr e-'"'^^\ (2.59) 

The latter relation is equivalent to 

^{lu; A) = $(cj; A) + Y r 5A$(a; X)da . (2.60) 

A J ^ 

Formula ( |2.60| ) is our key relation which together with (|2.56|) enables one to compute the 
physical spectral density ^{u) by using powerful heat kernel techniques. It is remarkable 
that (|2.60| ) reappears in the Euclidean theory (see section 3.3). 
Consider now a short t asymptotics of ( p.57| ) 



1 00 

Tre-*^^^') - ,,^,,<o-iy2 E «n(A)r. (2.61) 

n=0 



(47rt)(^-i)/2 



The coefficients an{X) are heat kernel coefficients for the operator ( p. 271) . They are local 
functionals on S^ and are even in A (the fiducial theory is U{1) invariant and the heat 
coefficients are even functions of charges). The gauge invariance also guarantees that the 
coefficients are polynomials in powers of the Maxwell stress tensor and its derivatives. 
In our case the role of the gauge field is played by the vector a^, and the corresponding 
Maxwell tensor Fij = aj^i — aj^i is related to the rotation. In general, 

[n/2] 

a„(A) = ^ X'"'a,^,n , (2.62) 

where a2m,n do not depend on A. The highest power of A in ( p.62| ) can be determined by 
dimensional analysis. Coefficients Oq and Oi in ( |2.61|) do not depend on A. 

The spectral density at high frequencies can be found from (|2.61|) by using ( |2.59| ). One 
can neglect for simplicity the mass gap and use inverse Laplace transform in (|2.59|) . It 
should be noted, however, that for operators with zero gap one has to take into account 



the presence of infrared singularities which come out in ( 2.57 ) at small 00. One of the pos- 
sibilities to avoid this problem is to use dimensional regularization and formally consider 
D as a complex parameter. 

It is instructive first to obtain the asymptotics for the fiducial spectral density at large 
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One can easily verify that for complex D substitution of ( p.63| ) in ( |2.57] ) yields ( |2.61| ). 



Another method how to define (|2.63|) for integer D is discussed in [Q. For ^{uj]X) 
relation ( p.59|) results in expansion 

a„(A) = a„(A) + —dxan+i{\)- (2.65) 



Finally, by taking into account (|3D, (|2:62|) , (|23|), (|2:65| ) one finds 



*<-'^(i^SF(#r;0"^"' '"" 



2n -P/™ 1?-1 



r [m 



Cn — 2^ (~1) ~;^ D-l\ '^2(m-n),m; (2.67) 



2 

We call Cn the spectral coefficients. It is remarkable that they are local functionals ex- 
pressed in terms of heat-kernel coefficients of elliptic operators. 

Now some comments are in order. Formula ( |2.67|) is valid in any dimension D. As for 



expansions ( p.63|) and ( p.66|) , they are finite in even dimension D = 2k. For D = 2k + 1 



terms in ( p.66|) with n > k are formally zero. This indicates that one should be more 
careful with infrared singularities in this case. The proper way of dealing with the infrared 



problem is to keep D complex in (|2.66| ) till the last stage of computations. Then for any 
D some relevant quantities determined with the help of $(0;) are finite except, possibly, 
a number of standard poles (see section 2.5). 

We return now to the question about asymptotics of the pseudo-trace K{t). To find 
its behavior at small t it is enough to substitute (|2.66|) in (|2.49|) . By neglecting the gap 

we get the expansion 

-1 00 

As was expected, (p^.68D holds both for odd and even D because infrared singularities play 



no role. We conjecture that the pseudo-trace expansion ( p.68|) has to be equally valid for 
continuous and discrete spectra. 

The coefficients in leading terms in ( p. 6^ ) can be immediately computed by using 
( |2.67|) . For instance, Cq = a^ and 

D-3 , , 

ci = ao,i H — a2,2- (2.69) 

Term 02,2 is determined by the gauge part of 02 (A), see ( p.62|) . 



02 2 = 

12 



^ / h'/'d''''xF^^F,„ (2.70) 



where Fij = aj_i — aij is a "Maxwell tensor" of the fiducial gauge field. Term ao,i coincides 
with Oi in expansion (|2.61|) for operator (|2.27|) with A = 0. 
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2.4 Examples 



There are simple examples where our results can be checked. One example is motivated 
by studying a quantum theory in the Einstein universe R^ x S^^^. Consider a Gibbs 
state in this space-time defined as a thermal equilibrium in the frame of reference which 
rigidly rotates with angular velocity Qq (with restriction Qq < 1/p where p is the radius 
of 5*^^^). The spectrum can be easily found for conformal scalar fields 



'-V^ + ^r)R)4> = 0, 



(2.71) 



where ^d = {D — 2)/(4(D — 1)) and R = D{D — l)/p^ is the scalar curvature. Suppose 
for simplicity that D = 3 and p = 1. The metric is 



ds^ = -dt^ + sin^ ed^^ + de\ 

In the rotating frame this metric can be written as 



„• 2 a 

ds^ = -B{dt + a^difi)^ -\ dif"^ + d6 



B 



where ip = ip — Vt^t and 



B = l-nl sin^ e, a^ = fio sin^ OB^K 



(2.72) 



(2.73) 



(2.74) 



The spectrum of frequencies in a non-rotating frame can be easily found, Un = n + 1/2, 
where n = 0, 1, .. is the angular momentum of a particle on 5*^. The frequencies in the 
rotating frame are ojni = ^n + ^o^ (see ( p..l8|) ), where |/| < n. Thus, the pseudo-trace 
K{t) is defined as 

oo n 

Kit) = Y.Y: e-^-*- (2.75) 

n=0 l=—n 



The first spectral coefficients in the asymptotic expansion ( p.68[ ) for ( p. 751 ) can be found 
by a direct computation^^ 



Co 



Air 



n'n ' 



Cl 



3 



ni 



(2.76) 



We show now that results ( |2.76| ) are reproduced by formula ( p.67| ) and cq, ci can be 
represented as integrals of local quantities. To apply (|2.67| ) we have to write the fiducial 
Hamiltonian in the form ( p.27| ) by using the conformal transformation ( p.26| ). Because 
the theory is conformally invariant we get|^ 



H\X) = -W^{V, + i\ai){V^ + iXa^) + \r+ ^F^^F.^, 

o oZ 



(2.77) 



^^For instance, one can define for this purpose the zeta function ^(z^) — J^ni ^ni " ^"^^ S^* "-o as the 
hmit 27r(i' — l/2)C(i/) at i' = 1/2 and ci as a hmit —Att{v + l/2)(^{v) aX v — —1/2. These hmits can be 
computed di rectl y from the series. 

^•^To get ( ^.77| ) o ne ha s to take into account that scalar curvature of 3D space with metric B^^ds^ 
where da^ is metric ( |2.72| ) is i? + ^F^'^Fik. 
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where ajrfx* = a^dip, Vj and R are covariant derivatives and the scalar curvature on the 
2- dimensional space B with the metric (see (p.73| )) 



dl' 



sin^^ 



1 



#" + -d9' 



52 "^ ' B' 

B and a^p are defined in (|2.74| ). For the operator (|2.77|) we get 

h}l^d^x = Vol B, 



Co 



B 



Cl 



¥/^d^ 



X 



B 



24^ 32 '' 



(2.78) 

(2.79) 
(2.80) 



where Fij = Fij is the Maxwell tensor for the fiducial potential Oj. By using (|2.74| ) one 
can check that ( p.79| ) and ( p.80| ) coincide with ( p. 761) . In the same way one can make the 
direct check of ( |2.67 ) for a conformal field in a rotating 4-dimensional Einstein universe. 



2.5 High temperature limit 

Let us consider a Gibbs state in a stationary space-time in the limit of high temperature^^ 
This limit is one of the most interesting and well studied regimes. When the parameter (3 
is small the dominant contribution to the free energy ( |1.6|) results from large frequencies 
u ~ j3~^^. In this situation one can go from summation in (|1.9|) to the integral 



/•oo , 

F^{p)=p-^ <^{uj)duj\n[l 



~/3w 



(2.81) 



where the spectral density ^{u) should be determined by its asymptotic expansion (|2.66| ) 
This gives 

1 



F^{(3) 



Er 



D-2n\ , . IP 



2n 



(2.82) 



Here ({x) is the Riemann zeta- function. High temperature asymptotics of this form have 



been obtained and studied on static space-times by Dowker and collaborators [|1^, [ ^ , 
(see also [^]). Equation ( 2.82 ) is the extension of these results to non-static case. 
The following remarks regarding ( p. 821) are in order. First, formula (|2.82|) is purely 



local. However, the method we used to get (|2.82|) does not take into account a non-local 
term proportional to /3~^, see fl^, |3^. This contribution comes for Bose fields from 



the region of small uj where the logarithm under the integral in ( |2.81|) results in a term 
ln{uj). Second, (|2.82|) is obtained with the help of dimensional regularization. When 
the parameter D coincides with the physical dimension one of the terms in (|2.82 ) has a 
simple pole. However if the system is in a finite box or there is a mass gap /x > this 

^'^Certainly, there must exist necessary physical conditions which ensure the thermal equilibrium. These 
conditions require special analysis in each particular case. We just assume that in the problems considered 
here they are satisfied. 
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singularity is an artifact of the method. Careful analysis of the infrared singularities has 
been carried out recently by Gusev and Zelnikov |J0| who have derived an approximate 
non-local expression for the one-loop free energy by using covariant perturbation theory 
of Ref. [^. Further discussion of this problem can be found in |[40|| . 
In four dimensions {D = 4) 



F^i(3) 



vr^ Co 



90 /?4 



24 /?2 



167r2 



ln(/5p)c2 



167r5/2 



n=3 



n-^)c(2n-3)cj^ 



2n-4 



(2.83) 



where p is a dimensional parameter related to the regularization and the pole term is 
omitted. When the frame does not rotate (|2.83| ) coincides with the result of Ref. [|14| . Of 
special interest are the leading terms in (|2.83|) 



F^{(3)^~Jd'xV^ 



90 24 



-R-V 



6 



-n' 



_C2_ 

16tt' 



HPp). (2.84) 



Here T is the local temperature T = l/[l3yB) and R is the scalar curvature of the 
background space-time. The second term in r.h.s of ( p. 84 ) depends on the local angular 
velocity Vt of the system with respect to a local Lorentz frame, see ( p.l2|) . This result 
follows from the form of the Ci coefficient ( p.69|) which can be rewritten in terms of the 
characteristics of the physical space-time as [Q 

2. 



Cl 



'—gdx 



B 



-R-v--n'^ 

6 3 



(2.85) 



As for C2 in ( p.84|) , we will show in section 3.5 that it does not depend on the rotation 
and has a standard covariant form quadratic in four- dimensional curvatures. 

The validity of (p.84|) can be checked for some simple models where the spectrum is 
known. The simplest case is again conformal fields in rotating Einstein universe. Such 



models have been recently investigated by several authors |^-|^5|]. The formula ( |2.84| ) 
is in agreement with these results. 



2.6 Thermal part of the stress-energy tensor 



The free energy F^{f3) of a Gibbs state on a stationary space-time is a functional of 
the background metric g^^p and the Killing vector ^. It is invariant with respect to time 
independent coordinate transformations because these transformations do not change the 
form of the metric, see (|2.6|)-(pl9|). Infinitesimal changes of coordinates can be written as 



^x^" = C^(x*), 



(2.86) 



where /i runs from 1 to -D and i from 1 to D — 1. The condition that coordinate changes 
are time independent can written as 



{c^cr = -(AO^ = 0, 



(2.87) 
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where C, = dtis the corresponding Kilhng field, and C^ is the Lie derivative along ^. Thus, 
( |2.87 ) is equivalent to the condition that the Killing field is left unchanged under trans- 



formations generated by (^^(x*). In a stationary space-time the group of transformations 
is D-dimensional and by the Noether theorem the invariance of F^{(3) implies existence 
of D Noether currents and corresponding charges. Like in a classical theory, the Noether 
current can be defined as T^i/i^^ where T^,y is the stress-energy tensor 

y„„ ^ 2 (^I^\ p.88) 

According to ( |2.87| ) the variations are to be taken at fixed ^ 0. Note that -F(/5) is a 
functional in a -D — 1 dimensional space, so g^^i, are functions oi D — 1 coordinates. One 
can use ( |2.84| ) to find T'^'^ in high-temperature limit As an illustration, we present the 
result for conformally coupled scalar field {V = R/6 in ( |2.84| )) 




+ — {2A,pA^ P - u.A^xw'' - u^A.y^w^ +u,A^^ '^ + u^A^^ •^) + O(lnT). (2.89) 

Here u^ is the velocity vector of the Killing frame where thermal equilibrium is defined, 
w^ is its acceleration and A^^, is the antisymmetric rotation tensor (|2.5|) . The invariance 
of F{j3) implies the standard covariant conservation law T^^^.^ = 0. Indeed, by considering 
F{P) as local functional ( p. 84 ) one can always write the change of F'^\P) under coordinate 



transformations as 

(^cF^(/3) = J d'x./^ [V,j^(C, 6^g) - T'^^Cm;.] , (2.90) 

where j^ is a local expression of ( and variations of the metric. Because the background 
theory is stationary the divergence of j^ reduces to three-dimensional divergence which 
results in a surface term. For stationary transformations ( |2.8tj| ) the expression (2.9U) 



vanishes which is possible only when T^^" is divergence free. This justifies its interpretation 
as a stress-energy tensor. 

By taking into account that u^A^".^ = —A^^yA^^ one can check that T/^ = 0. Hence 
is interpreted as a stress-energy tensor of a massless conformally invariant matter 



with temperature T and velocity u^. The first term of ( p^.89D coincides with the stress- 
tensor of an ideal gas at temperature T. 

It should be noted that (|2.88D represents the average value of the thermal part of a 



quantum stress tensor-energy tensor in the Gibbs state. It does not coincide, however, 
with the total stress-energy tensor because ( |2.88[ ) does not include the vacuum partQ 



^^In practice one has to keep fixed contravariant components ^'^ of the Killing field. For instance, 
variation of i? = — ^^ then looks as {5B)^ = "^^^"^giiv- 

^^Let us recall that the contribution of zero-point energies in F'-^ {j3) has been eliminated from the very 
beginning by normal ordering, see section 1.1. 
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Because the conformal anomaly is related to renormalization only of the vacuum energy 
this fact explains why ( p.89|) is traceless. 



High-temperature asymptotics like ( |2.89| ) may represent an interest in studying quan- 



tum effects near rotating black holes because the local temperature T becomes large near 
the black hole horizon. To realize such a state for a Kerr-Newman black hole one has to 
place it in a cavity which co-rotates rigidly with the black hole. Hence, for Kerr-Newman 



black hole ( p.89|) requires to take into account effects of the cavity. The situation is simpler 



for rotating black holes with asymptotic anti-de Sitter geometry (Kerr-AdS black holes 



2%). They are solutions in gravity theory with negative cosmological constant. If the 
angular velocity of Kerr-AdS black hole is less than inverse AdS radius the Killing vector 
field ^ is globally time-like outside the horizon. The cavity is not necessary in this case 
and (|2.89|) may be a correct approximation to the thermal part of stress-energy tensor, 
at least near the horizon. 



3 Euclidean theory 

3.1 Formulation of the problem 

Let us consider a stationary space-time A4 with metric defined by (|2.6| )- (|2.9| ). Suppose 
that the metric depends on a number of parameters J denoted by a collective symbol J. 
These parameters will play the same role as the angular velocity in the metric ( |1.17|) . We 
assume that there is a complexification of Ai (in the sense described in section 1.2) which 
has a real Euclidean section Me with metric 

ds% = B{dT + dkdx^f + hkjdx^dx^ . (3.1) 

r is a periodic coordinate with period /3. We also assume that (|3.1[) can be obtained from 



by replacement t = —it and analytical continuation of the components to imaginary 
values of J 

hkj{x; J) = hkj{x;iJ) , 

dk{x; J) = iak{x; iJ) , B{x; J) = B{x; iJ) . (3.2) 

An example of this procedure has been discussed in section 1.2. The Euclidean effective 
action W^{f5, J) for free fields or in one- loop approximation is a determinant of the 
corresponding Euclidean operators Le, see ( |1.1(]| ). In what follows we consider the action 
defined with the help of ( function regularization P6|],[B^ as 



WiP, J) = -\ lim |;C(^l/5, J). (3.3) 

C(H/3,^) = E[^'A]"^ (3.4) 
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where A are eigen- values of Le, £> is a dimensional parameter related to regularization. 
The series in ( |3.14 ) converge for Re v > D/2 — 1. To get ( function for other values of 



z/ one considers analytical continuation. It can be shown that ({I'lP) is a meromorphic 
function of z/ with simple poles on the real axis. In particular, this function is regular 
near z/ = so that one can use ( p.3|) . 

To go to the physical space-time one has to analytically continue (|3.3| ) back to some 
real values of J (these values are to be in an interval corresponding to the rotation with 
the velocity smaller than the velocity of light). This results in a new functional 

f3F^{(3) = W^{(3,-zJ). (3.5) 

We call F^{j3) the Euclidean free energy. Our purpose now is to describe the relation 
between F^{[3) and canonical free energy F^{[3), Eq. ( |1.9| ) 

As earlier we consider the scalar model ( ^2.181 ). Then (^-function ( |3.4| ) is determined by 
the eigen- values of the Euclidean wave operator on A^^; 

(-V^V^ + V)<Pk = A<j)A. (3.6) 

The Kaluza-Klein method can be applied to Euclidean theory as follows. Note that 
because of the isometry, 0a are eigen-vectors of the operator idr 

0A(r,a;)=e*'^'>A(a:), (3.7) 

where ai = {2ttI)/i3, I = 0, ±1, ±2, .... For these modes ( |3.6D takes the form 

{-9Ei^^. + icria^){^u + i^ia,) + V)(Pa = Hk- (3.8) 

Here the metric and the connections correspond to a fiducial static Euclidean space-time 

Me 

ds% = {gE)fiudx^dx'' = Bdr^ + hjkdx^dx^. (3.9) 

Thus, the Euclidean problem on a stationary background can be reformulated as a theory 
on a static background in the presence of a fiducial gauge connection. The distinction 
from the Lorentzian theory is that the fiducial charges ai (Matsubara frequencies) are 
quantized because the Euclidean time coordinate r is compact. 

Without loss of generality we again suppose that _B = 1 in ( p.9| ) . One can always use 
conformal transformation to bring the metric to this form. The result of this transforma- 
tion is an anomalous addition to the Euclidean action. This addition is proportional to 
/3 and changes only the vacuum energy. We will discuss this question in section 3.4. If 
B = 1 the eigen-value problem can be written as 

{af + H\ai))(j)A = A(j)A, (3.10) 

H\ai) = -h^\Vj + iaid^){Vk + i<yidk) + V. (3.11) 
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The operators H{ai) are analogous to the Lorentzian Hamiltonian H{X). H'^{ai) are 
positive-definite elhptic operators. 

Let us consider first systems on a compact space. The spectrum of n{(Ti) is discrete 
and bounded from below and we suppose it is strictly positive, for simplicity. Denote the 
corresponding eigen- values by uj'^{ai). Given solutions to the problem 

H\a{)t^\x^)=u\a{)cl^t\^'). (3.12) 



the solution to ( |3.6| ) can be represented as 

0A(r,x) = e^'^'>(r')(xO, A = a^ + uj\ai). (3.13) 

Then (^-function (|3.4| ) is 



ai u} 



where the sum is taken over all eigen- values Lo'^{(yi) of H'^{ai) (some eigen- values may 
coincide). The Euclidean effective action is defined with the help of (|3.14|) by ( |3.3|) . 

3.2 From Euclidean to Lorentzian theory 

By analogy with the Lorentzian theory, where fiducial charges are continuous, let us 
consider a one-parameter family of operators -ff^(A) where the parameter A is real. For 
the considered problem H'^{X) are positive-definite elliptic operators with eigen- values 
u'iX). 

We make the following additional assumptions: 

1) The Wick rotation to Lorentzian theory is determined by analytical continuation 
of a set of parameters J in such a way that under this continuation eigen- values of H^{X) 
transform into eigen- values of the Lorentzian operators H'^{X), i.e. 

u\iX\-iJ)=uj\X\J). (3.15) 

Equation ( p.l5| ) should be considered simply as a condition on lj"^ as functions of A and 
J. 

2) As in the Lorentzian theory, see ( p.32|) , we define functions x(cij. A) = cij^(A) + A^. 



For X = ai they coincide with eigen-values A of the Euclidean operartor Le- We assume 
that x{u!, z)~^ can be analytically continued in z to the upper complex plane where they 
are meromorphic functions with simple poles. Note that according to (|3.15|) , x(u;. A) 



transform to — x(to'. A), functions defined in ( p. 321) . Thus, after the Wick rotation the 
poles lie on axes Re z = and coincide with physical energies. The requirement that 
the poles are simple is equivalent to condition ■x'ioj^iuj) ^ 0, where a prime denotes the 
derivative with respect to A. This is in accord with condition that physical states have 



non-zero norm, see ( p.44|) 
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For simple models one can check that our second assumption about x^u, A) does hold 
true. For instance, for a field in a rotating Einstein universe (discussed in section 2.4) the 
eigen values of the Euclidean operator Le = — V^ + -R/8 are A = {ak + ICIqY + u^ where 
CIq is the Euclidean angular velocity. Hence, x{uj, A) = (A + ICIq)"^ + u^ and under change 
CIq to —iQo one gets x(cc;,'iA) = —(A — IQq)^ + uj^. They are the functions whose zeros 
A = ±{ujn + /fio) coincide with physical energies. Also, x'{uj, ioj) = 2iujn ^ 0. 



We now follow |48] and rewrite the (^-function by using the Cauchy theorem as 



-2v 



C(H 



-EE 



27ri ^^Jcz-ai 



dz 



{z' + iJo\z))- 



(3.16) 



The contour C consists of two parallel lines, C+ and C_, in the complex plane. C+ goes 
from {ie + oo) to {ie — oo) and C_ goes from (— ie — cxd) to (— ie + cxd). Here e is a small 
positive parameter such that e < fi. We consider (|3.16|) at Re u > D/2 — 1 and assume 
that D >2. Summation over ai in ( |3.16|) can be performed 



cW)-^E//-ot(^^ 



{z' + u;\z)y 



(3.17) 



Note that spectrum (i^^(A) has to be symmetric with respect to change A to —A. Hence, 
integrations over C+ and C_ in ( |3.17| ) coincide and one can write ( p. 17] ) as twice the 
integral over C+. To proceed let us use the identity 




pdz V 



Jf3z 



which enables one to write 



C(H/5) = /?Co(^) + Ct(H 



(3.18) 
Co(z/) = ^E/ {x^ + ^\x)r^dx, (3.19) 

Ct(H/5) = ^-E/ dz^ln{l-e^^^){z' + u\z))-r (3-20) 

Ct(^|/?) represents the purely thermal part which vanishes at zero temperature because 
of small positive imaginary part of z in e*^^. This means that the only quantity which 
can be responsible for the vacuum energy is Co('^)- Let us now deform C+ in ( p. 191) so as 
to make the integrand exponentially small at large z due to the factor e*'^^. After that we 
can integrate by parts to get 



Ct(^| 



Q 



-2v 



TTl 



Y, [ dz\n{l-e'^' 



XK^.z) 



(3.21) 



IC+ ^ ' {z^ ^ uj\z)y^^' 

We can now represent the Euclidean action in the following form (see ( |3.3| ), ( |3.18| )) 

W^^(/5)=/3(F(/5) + Eo), (3.22) 
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F(/5) = -llimAcTW), (3.23) 

where Ct(i^|/5) and Coi.^-') are given by ( p.l9| ) and (|3.21|) , respectively. To compute F{[3) 
we note that Cr(i^|/5) has a form Cl'^l/?) = ^fi^lP)- To find /(i^|/5) at z/ = we add 
to C+ a large semicircle lying in the upper half of the complex plane to make a closed 
contour. Because of the exponent e*^^ in the logarithm in ( p.21| ) the integration over the 



semicircle vanishes when its radius goes to infinity. Our assumption guarantees that the 
function /(z^|/5) is finite at z/ = and by the Cauchy theorem its value is determined by 
the residues at z = iu. We get 

W) = ^Eln(l-e-'^^), (3.25) 

where z are corresponding solutions to x('^; ^-2) = 0. Finally, the first assumption guar- 
antees that after the Wick rotation ( p.25| ) coincides with the canonical free energy 

F''{f3) = W^n{l-e-^-), (3.26) 

where the sum is taken over all non-negative solutions of the equation uj'^{z) = z^ (see 
( p.31| ), ( p.32| )). Thus, according to ( p.5|) , ( p.22|) the canonical and Euchdean free energies 



are related by Eq. ( |1.16 ) where Eq coincides with ^0 after the Wick rotation. We will 



see in section 3.4 that Eq does correspond to the vacuum energy. 
3.3 Continuous spectrum 

It is instructive to discuss what happens in case when the spectrum is continuous. The 
Euclidean theory offers another look at the origin of key relations ( p.56| ), ( |2.60|) which 
determine the physical spectral density (^{u). 



The Euclidean action is defined by (|3.3|) where the C function now has the form 

CW) = g'"' dooY.^^^^di^f + ^'r"- (3.27) 

We assume that n{(Ti) have a positive mass gap /i which does not depend on Matsubara 
frequencies o"/. Similarly to the discrete case, the (^-function can be rewritten by using the 
Cauchy theorem in ( p.l8| ) where 



Co{iy) = - / duj <l>{uj;x){x^ + uj^)-''dx, (3.28) 

vr Jfi JO 

(^{j,\f]) = ^— duj dz— In (l-e'^')<!>{ij;z){z^ + ijY\ (3.29) 

TCI J fi Jc+ dz ^ ' 
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and C+ is defined as earlier. Then the Euchdean action is represented in the form ( |3.22| ) 
where F{i3) and ^o are defined with the help of ( |3.28| ) and (|3.29| ), respectively. 

To compute Ct{^\P) we proceed as before and deform C+ in ( |3.29|) so as to make the 
integrand exponentially small at large z due to the factor e*^^ |^. Integration by parts 
over z then gives 



Q' 



-2u 



oo 



(t(u\P) = - / duj dz\n{l-e"^' 

T^l Jfi Jc+ ^ 



2uz^{uj;z) d;,^{uj;z) 



(3.30) 



One can also integrate by parts over u the second term in square brackets^ 

CtW) = '^^ r du; f dz ^^\ \n{l-e^^^)n^;z), (3.31) 

m JfM Jc+ [z^ + uj'^y^^ ^ ' 

^(w; z) = ^{u; z)-- r ^^{a; z)da. (3.32) 

z Jfi dz 

Function '^{uj]z) corresponds to auxihary density \Ef(ti;,A), see (|2.60|) . In the Lorentzian 
theory \l'(c(j, A) appears as result of analysis of the inner products. It is remarkable that 
in the Euclidean theory there appears an analogous function ^(tu; 2;), although it happens 
on a completely different footing. 

To compute F{(3) one can use the fact that Ct(i^|/5) now has a form C(^|/^) = ^f{^\P)y 
see (|3.31| ). We follow [|^ and assume that for any u the density ^(cu, z) can be analytically 



continued to complex z and is an entire function of z in the upper half of the complex 
plane. One can then proceed as in section 3.2 and compute f{i^\P) at z/ = by using the 
Cauchy theorem. The result looks as follows: 

^ 1 roc V / \ 

F{P) = - / du^u) In (1 - e-^^j , (3.33) 



^(00) = '^(uj\z] 



(3.34) 



Similarity between <l>(u;) and physical density ^{u) defined by (|2.56|) , ( |2.60|) is obvious. 



Consider now the Wick rotation from Euclidean /Ae to Lorentzian space-time A4 deter- 
mined by analytical continuation of a set of the parameters J of the metric. Our second 
assumption is that ^(cu; A), as a function of parameters J, can be analytically continued 
to complex J in such a way that the following equality holds: 

^(cj;iA| -iJ) = ^(cj;A|J) (3.35) 



for J and A real. Condition ( |3.35| ) is an analog of ( |3.15| ). Given equations ( p.33| )-( p.35| 



one concludes that ^{u) and F{P) coincide with the physical density of levels and the 
free energy, respectively, 

^uj\-iJ) = ^{uj\J), (3.36) 

^^Wc conjecture that $(tj; z) itself cannot increase at large z and presence of the factor e^^^ is enough 



to ensure convergence of ( 3.29 ) . 

^^The boundary terms vanish at Re ;/ > D/2 — 1 because dz'^{uj; z) ~ lli^^"^ at large w. 
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F{P\-iJ)=F{P\J). (3.37) 

A rigorous proof of ( p.35| ) is problematic because the explicit form of $(u;; A) is not 



known in general. One can demonstrate [Q, however, that this assumption is true at 
least asymptotically in the limit of high frequencies when the spectral density becomes 
a local functional of the background geometry and the Wick rotation procedure can be 
easily applied. 

3.4 Vacuum energy 

Let us return to the vacuum part ^o of the Euclidean action, see (|3.24|). According to 



(H), 



io = -^lim-^Co(^^), (3.38) 

Uu) = — dxC Ug\H\x) + x')) , (3.39) 

ng Jo ^ ' 

where (^{v\0^ is the generalized ^-function of an operator O. Let us show that after the 

Wick rotation Eq agrees with the standard definition of the vacuum energy. For this 



purpose we represent (|3.19| ) in the form 

Co(^) = -^ E / ^^(^' + '^'(^))-^e-^ (3.40) 

27r V-^c+ 

where a small positive parameter e is introduced to regularize the integral in the limit of 
vanishing v. By integrating in ( p.40| ) by parts over z^ and neglecting terms linear in e, we 

get 

The regularization enables us to replace C+ by a closed contour in the upper half of 
complex plane and use the Cauchy theorem 

^o = ^E^e-^ (3.42) 

where z are solutions to xi'^^ ^^) = 0- After the Wick rotation Eq does correspond to the 
vacuum energy determined with cutoff 1/e in the range of high frequencies. 

A comment regarding relations (|3.22|) and (|3.42|) is in order. These relations are valid 



when the Killing field ^ has a unit norm, i.e., the space-time is " ultrastationary" with 
metric 

ds'^ = {dr + dkdx^Y + hkjdx'^dx^, (3.43) 

where r is periodic with period (3. One can always use the conformal rescaling g^^iy = g^v/ B 
to reduce the metric g^^, see (|3.1|), to ultrastationary form ( p.43| ). The (renormalized) 



Euclidean action VT'^ff?] on ( p.l|) is related to the action VT'^f^'] on (|3.43| ) as 

W'^[g] = W''[g]+p/\[gi (3.44) 
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where /3A[5f] appears as a result of the confomial anomalyQ. It can be shown that A[g] 
is a local functional [^,|3^. According to ( |3.43|) , the vacuum energy computed for an 



arbitrary stationary space-time is Eo[g] + A[g], where EqIq] is defined by ( p.38|) , (|3.39| ) on 
on (lOSI) 

Equations (|3.38| ), ( p.39| ) hold for discrete and continuous spectra. They enable one 



to compute the vacuum energy of fields in arbitrary stationary space-times by using the 
(■-function method and can be used in different applications. 

3.5 Kaluza-Klein reduction of heat-kernel coefficients 

Let us discuss now the functions ({ulf]), Col'^) and relation (|3.18| ) (which holds on ul- 
trastationary space-times) in more detail. We begin with ({I'lP) = C {'^Iq'^L) which is a 
generalized ( function of the Euclidean operator Le = — V^ + V, see ( p. 61) and (|3.14|) . 
According to the general theory, ({I'llS) is a meromorphic function which has simple poles 
on the real axis of u. In the theory with D dimensions the part C^^\i^\P) which includes 



the poles of Cl'^l/^) can be written as (see, e.g., [|17| ) 



2 °° A 

where An are the coefficients related to the asymptotic expansion of the heat kerneQ 

Relation ( p.l8|) has an important consequence: the poles of ((^^1/^) coincide with poles of 
j3(o{i'). This happens simply because coefficients An in ( |3.45| ) depend on parameter /? 
hnearly. To investigate the poles of Co('^) we rewrite ( |3.39|) as 

(o{u) = — — / dx dt r-^Tr e-*(^' (")+"'). (3.47) 

7rr(z/) JO Jo 

Our assumption is that the poles of (q are related to the behaviour of the integral at small 
t. The argument is that in this limit the trace is a local functional. By following |l^ we 
define the pole part Co A'^) cif Co('^) by (|3.47|) with integration over t taken in the interval 



(0,1). One can then replace the trace of the heat kernel of H'^{x) by its asymptotic 
expansion 

-1 oo 

T--*"^^^^ - .^^,.n-,y. E UXr. (3.48) 

V / n=0 

^^The anomalous term is related to breaking of conformal invariance in the process of renormalization 
of the effective action. The anomaly can be also attributed to non invariance of the integration measure 
in the Euclidean functional integral | p9[ . 

^°In what follows we put p = 1 for simplicity. We will also assume that space-time has no boundaries 



and hence n in (3.46) is an integer. 
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and get 



d^V) 



oo /•! 

dx 



dt r-^e"*^' 



7rT(i^)Jo ^^Jo "^^ " (47rt)(^-i)/\^o 



E«n(a;)t". 



(3.49) 



The integral exists at Re z/ > (-D — l)/2. To proceed we note that the heat kernel 
coefficients are polynomials analogous to (p.62|) 



[n/2] 

,(A) = ^ A2-a2, 

m=0 



'm,m 



(3.50) 



where a2m,n are some coefficients. It enables one to integrate ( p. 491 ) over x and then over 



d^V) 



[n/2] 



EE 



r(m + l/2)a2m,^ 



(47r)^/2r(z/) „t^o^o v^(2z/ + 2(n - m) - D) ' 
The latter equation can be rewritten also as 



d^V) 



(47r)^/2r(z/; 



E 

n=0 



1 



2u + 2n- D 



^ r(m-r2 + l/2) , 

/ ^ / — '^2(m— n),m- 



TT 



(3.51) 



(3.52) 



By comparing poles of ( 3.45 ) and (|3.52 ) one concludes that 

^n = ^ E r("^ -n + l/2)d2(m-n),m- 



(3.53) 



For D odd this conclusion is true for all n while for D even only for < n < D/2 — 1 
(terms with other n do not result in poles). It is clear, however, that ( p.53|) should be 
universal for all D and all n because dimensionality does not appear in it explicitly. 

Equation (|3.53 ) can be considered as a Kaluza-Klein reduction formula for heat kernel 
coefficients of an operator on a stationary D dimensional manifold. It has one interesting 
and important consequence at n = D/2 for even D . Namely, quantity l3~^A£,/2 transforms 



under the Wick rotation to the spectral coefficient C15/2 defined by ( |2.67| ). This can be 
easily seen with the help of ( |3.53|) if we take into account the relation between Euclidean 
and Lorentzian coefficients 

a„(— iJ, iA) = a„(J, A). (3.54) 

Equations (g3§, ( ^TU^ and (^) imply that 



d2m,n{-iJ) = (-l)"'a2m,n(^)- 



(3.55) 



The fact that cd/2 is related to l3~^A£,/2 is remarkable for two reasons. First, because 
Ad/2 = (4vr)'^/^C(0|/?) and C(0|/3) determines the anomalous scaling of the Euchdean 
effective action [^. Second, because this fact implies that cd/2 is a local covariant 
functional of D-dimensional curvatures. One can, thus, conclude that the coefficient C2 
which appears in high temperature asymptotic ( 2.84 ) in four dimensions is a covariant 
functional quadratic in curvatures, the same property it has in static space-times. 
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The argument we used to derive ( |3.53| ) is not quite rigorous. Relation ( |3.53| ), however, 
passes a non-trivial check. Consider (|3.53|) for ra 7^ on a D-dimensional stationary 
Euclidean background ( p.43|) where r is periodic with period /?. We denote this space by 
J^E- The D — 1- dimensional space B corresponding to operators H^{x) is determined by 
the metric 

df = hkjdx^dxK (3.56) 

We denote the Riemann tensors on A^^; and B as R^^xp and R^uXp, respectively. Their 



relation is discussed in Ref. 48 



Consider (|3.53|) for n = 1. It can be shown |^ that in any dimension 
Ai = pJ Vhd^'-'x i^-R-v'^=P j Vhd^'-'x (^^R - Y^F'^^^o.p - V^ , (3.57) 

where the last line holds on ultrastationary space-times. On the other hand, according to 



( p33D 



Ai= p ( ao,i + -02,2 

By using heat kernel asymptotics in external gauge field one finds 

-1 



(3.58) 



Vhd 



D-l 



ao,i = / V na" ^x 1 -n ~ v j , (3.59) 

02,2 = - / Vkd'^-^xj^F'^^Fap. (3.60) 

Hence, the right hand sides of (|3.57| ) and ( |3.58|) do coincide. Equation ( p.53|) can be also 
checked for the coefficient A2, which on a manifold without boundary looks as 

A2 = pjVhd^-'x Q^Pe - ^^KuR'" + ^^K^x.R'"'"' - \rv + 1^^) (3.61) 

The relation ( p.53|) for n = 2 is 



M= I3\ do,2 + 2^2,3 + 7^4,4 



(3.62) 



The coefficients ao,2, 02,3, 04,4 can be found by using results of Refs. |5^, [^, i.e 



0-0,2 



/ ^/hd°-'x (—R^ - —R.^R'"' + —RauXpR""^' - - W + ^V") , (3.63) 
i V72 180 ^ 180 ^^ 6 2 / ' ^ ^ 



02,3 



\fhd 



D-l 



12 M^^ ^ go ^''P 



A 



||A - -{FP''F,,)R 



Qj4. 4 



180 



.^Ai^^p^ P 



\ii^ fj,p 



R PPM p V 

90 ^'^ ' 



288 



1 
3"6b 



{FP^F^^f + F^^F^PFy^^F.p 



(3.64) 
(3.65) 



Here the symbol || corresponds to covariant differentiation with respect to the metric h^^y. 
By using (g:62D-(F65|) one can check that ^M^ holds |8|. 



Formulas (p.53|) are also interesting in applications to genuine Kaluza-Klein theories. 
In this case the Euclidean space M.E has to be identified with a higher dimensional space, 
the Euclidean time with an extra coordinate and /? with a compactification radius. 
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4 Generalization and external gauge fields 

As we have seen the Euchdean formulation of statistical mechanics offers an alternative 
derivation of relations ( p.56| ), ( |2.6CI| ) between physical and fiducial spectral densities. The 



advantage of this derivation is that actually it does not depend on the form of the Eu- 
clidean wave operator Le- Hence, one can conjecture that ( p.60| ) is more universal and 
may be true (with some additional assumptions) when Le is a higher-order elliptic op- 
erator. If this were the case such results of Section 2 as spectral asymptotics could be 
extended to non-linear spectral problems of a generic form. We will not analyze this 
possibility here but mention another problem where ( |2.56|) and (|2.60|) are satisfied. This 
problem is interesting for physical applications. 

Consider a charged scalar field in Minkowski space-time which interacts with a static 
electric potential ip{x'^). The equation of motion is 

r'^r"^ + m^) = 0, (4.1) 

where P^ = d^ — igA^, Afj_dx^ = ipdt and g is the charge of the field. For single-particle 
excitations (/)^(t,a;*) = e~**^0^(x*) is reduced to the non-linear spectral problem 

{uj''~H^ico))<PUx') = 0, (4.2) 

H\uj) = -df -g^ip"^ -2ujgip + m^. (4.3) 

There must be some physical restrictions on the potential if so as to avoid complex energies 
u in (^4.2| ). Complex u correspond to an instability of the system (creation of particle- 
anti-particle pairs, for example). To exclude such effects we assume that electric field is 
sufficiently weak and vanishes at spatial infinity. Note that equation (|4.1| ) preserves its 
form under the gauge transformations 0' = e~*^'^*0, ip' = (f + a, where a is a constant. 
These transformations shift the spectrum of u in ([4.2| ) to uj + ga. We will eliminate this 
arbitrariness by requiring that if is zero at infinity. 

One can now proceed as before and formulate a corresponding fiducial problem 

The Klein-Gordon inner product for physical and fiducial fields are determined, respec- 
tively, by the currents 

jMu 02) = -^{<PlV,<j)2 - {VM*<P2), (4.5) 

JMi, 02) = -^(015^02 - 5^0t02). (4.6) 

It should be noted that for charged fields one has to consider two Hamiltonians: one is 
if (a;) defined in ( |4.3| ) and the other is H{—uj), for u; > 0. The reason is that the system 
contains particles and antiparticles which have different charges and interact with the 
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electric field in different ways. The particles and antiparticles are described by H{uj) and 
H{—u), respectively. The total spectral density ^{uj) is a sum of the spectral densities 
of particles, <l>_|_(u;), and antiparticles ^^{u). The fiducial problem for antiparticles has 
to be formulated as ( [4.4| ) with A replaced by —A. Suppose that ^±{uj;X) are fiducial 
spectral densities for particles and antiparticles and \l'±(ct;; A) are corresponding auxiliary 
densities defined by ( p.60| ). Then by using ( |4.5| ), (f4.6|) one can prove relation ( |2.56|) , i.e. 
that $±(u;) = '^±{uj;uj). The asymptotic expansion for the total density $+(ct;) + ^-{u) 
is given by ( p.66|) , ( p.67|) . The spectral coefficients are c„ = c+^„ + c_^„ where c±^„ are 
determined by the heat kernel asymptotic expansion of operators H'^{X) and H'^{—X). 
These results can be used to show that an electric field yields a correction T'^ip'^ to the 
high-temperature asymptotics ( |2.84|) . 



One can come to the same result by considering Euclidean theory. Contributions 
from particles and antiparticles in this case appear in Euclidean ^-function, (|3.16|) , as 
integrations along contours C+ and C_ lying in the upper and lower parts of the complex 
plane, respectively. For charged fields these integrations do not coincide. 

5 Quantum fields near black holes 

5.1 Motivations 

As we have already discussed, quantum fields in thermal equilibrium with a black hole 
near its horizon appear, to a static observer, as a system at high temperatures. The free 
energy of such fields is described by high-temperature asymptotic in a form like ( [2.841 ). 



The local temperature T becomes infinite at the horizon. To avoid this divergence the 
integration in ( ^.84| ) has to be stopped at some small (proper) distance / near the horizon. 
It is easy to estimate that the entropy of the gas is of order of A/P where A is the surface 
area of the horizon. As was first pointed out by 't Hooft |S^, if / is associated with the 
Planck length Ipi the entropy of the gas is of the same order as the entropy of a black 
hole S^^ = A/ AG (here G = Ipi is the Newton constant). 

The entropy of black holes S^^ was introduced by Bekenstein |^3[ and Hawking pi 



who used the fact that black holes possess properties similar to the properties of ther- 



modynamical systems |^. In a classical theory, however, a Schwarzschild black hole is 
nothing but an empty space with a strong gravitational field. It still remains a fundamen- 
tal problem how to identify microscopical degrees of freedom of a black hole responsible 
for S^^ . The observation by 't Hooft is important because it offers one of the possible 
explanations of the Bekenstein-Hawking entropy. 

We do not give here a detailed analysis of this idea and do not include a (still growing) 



list of related references. This can be found in a review paper [Q . In remaining sections 
we discuss a single problem, i.e., the correspondence between canonical, F'^ , and Eu- 
clidean, F^, free energies in external regions of black holes. In particular we discuss the 
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duality property between infrared thermal divergencies of F^ and ultraviolet divergencies 
of F^ caused by conical singularities. This is an interesting feature which is important in 
finding connection between S^^ and the entropy of a thermal atmosphere around a black 
hole. 



Note that more detailed analysis of this problem is presented in Refs. [0,|^. Our 
aim here is to sketch the main results and pay attention to features related to rotation of 
a black hole. 

5.2 Killing horizons 

We begin with definitions. Consider a space-time M. which possesses a one-parameter 
group of isometrics parametrized by t and generated by a Killing vector field ^^. Suppose 
the isometrics leave fixed each point of a smooth space-like two-surface S. The surface 



S is called a fixed point set of ^^. It can be shown [^ that existence of S implies the 
existence of a Killing horizon on M.. The Killing horizon consists of two hypersurfaces 
spanned by null geodesies and the Killing field is tangent to these geodesies. The two 
hypersurfaces orthogonally intersect at S, and that is why E is a bifurcation surface for 
orbits generated by i^. The structure of these orbits near S is similar to orbits of Lorentz 
boosts in Minkowski space-time. 

The event horizon of black hole solutions in Einstein theory is the Killing horizonp^ 
The Killing horizon divides the space-time into four regions: two regions JF and V where 
^ is space-like, and two causally independent regions £ and TZ where ^ is time-like. In 
case of black hole space-times one of such regions, say 7^, describes the exterior region of 
a black hole. In such a region near E the metric can be written in the form 

ds^ = -B{dt + aadx'^Y + dp^ + aapdx^dx^, (5.1) 

where x" is a set of D — 2 coordinates and S is the surface p = 0. The coordinate p is the 
distance measured along a geodesic normal to S between a point with coordinates p, x" 
and a point on S with coordinates x". The components of ( ^.1|) do not depend on t. One 
can show that B ~ K^p^ at small p. Here n is the surface gravity, a constant which can 
be also defined in a covariant way in terms of the Killing field 






(5.2) 



The Euclidean space J^^ corresponding to A^ is obtained by the Wick rotation as 

we described earlier. The analog of S on A^^ is a hypersurface which is a fixed-point 

set of the Euclidean Killing vector field C,^ . The Euclidean metric can be written in a 

form similar to (|5.1|) where t is replaced by —ir. One can easily see that there are conical 

^^Note, however, that we do not require in this Section that the metric obeys the Einstein equations. 
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singularities on E if the period of Euclidean time r is arbitrary. A^^ is regular if the 
period is 27r//€. 

The geometry near Lorentzian or Euclidean horizons can be characterized by geometric 
invariants obtained by projecting components of the Riemann tensor and its derivatives 
on directions orthogonal to S. The corresponding projector is defined in terms of two 
vectors / and p normal to S as P'^^ = I'^l" =F P/^Pi/, where /^ = 1, p^ = ^1, {I ■ p) = 0. 
Signes — and + correspond to Lorentzian and Euclidean theoriues respectively. 

5.3 Canonical free energy 

In the presence of a Killing horizon, the spectrum of single-particle energies u has a 
number of specific features. To better understand them let us first return to the high- 
temperature asymptotic ( |2.84| ). The only difference between stationary and static space- 



times in this formula is the presence of the term T^Q"^ which appears when the Killing 
frame rotates with the local angular velocity Q. By using the definition of Q (see ( p.l2|) ) 
and the asymptotic form B ~ K^p^ in ( |5.1| ) it can be easily shown that Q^ ~ p"^ near the 
horizon and, hence, the term T^Q'^ is finite at p = 0. This means that the contribution of 
rotation in ( |2.84| ) can be neglected as compared to other terms (proportional to T^ and 
T^). Hence, fiducial potential a„ in ( |2.81| ) can be ignored and stationary space-time Ai, 
Eqs. (pl6|)-(pl9|), can be replaced by a fiducial static space-time /A, Eq. (|2.23|). Really, 



near the horizon one has effectively a static problem [Q . 

The last statement can be also formulated in geometrical terms. Consider three in- 
variants P^'^P^f'R^xup, P^^R^JLU1 and R computed at S, where R^^Xp, R^iu and R are 
D-dimensional Riemann, Ricci and scalar curvatures, respectively. It can be shown [Q 
that these quantities computed on JH coincide with corresponding quantities computed 
on Ai. Let us also note that surface gravities k of the Killing horizons as well as the 
bifurcation surfaces on A4 and Ai are identical. 

To understand what happens near the horizon let us neglect for a moment the curva- 
tures near E and approximate the metric by the metric on the Rindler space 

ds'^ = -n^p^df^ + dp^ + dzl + ... + dz]^^r^. (5.3) 

As was explained in section 2.2. the single-particle Hamiltonian can be obtained with 
the help of a conformal transformation which turns ( |5.3| ) into an ultrastatic space. From 
( p:27D , ( p30D one finds 

H"^ = -ViV* - aln^ + kVV. (5.4) 

Here Vi are covariant derivatives on the space B with the metric 

df = K'"^ p-"^ {dp^ + dzl + ••• + dzD-2)- (5-5) 

i3 is a hyperbolic (Lobachevsky) space with constant negative curvature R = —k,'^{D — 
1){D — 2). Although to leading approximation the potential V in ( p.4|) can be neglected, 
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H^ has an effective tachyonic mass —0]-,^? where ar, = {D — 2)/2. In its turn, — V*Vj 
on B has a mass gap which appears due to constant curvature. What happens is that the 



tachyonic mass exactly cancels the mass gap (see |5^,[^ for details). Thus, the spectrum 
of H^ is continuous and without mass gap. 

Because ( ^.5|) is a maximally symmetric space the spectral density of H^ grows as the 
volume of B. The volume divergence is of the infrared type and one can show that it 
appears in high-temperature asymptotics (|2.84|) as divergence of the T^ term. 

In previous sections our approach to volume divergences was to cut off integrations 
over the space to have a system in a box. In case of black holes it means that the region 
of physical space-time whose proper distance p from the horizon is smaller than some 
length e should not be considered [^ . The drawback of this procedure is that it makes 
the space-time incomplete. Thus, it makes sense to consider other regularizations which 



allow one to work on an entire external region of a black hole. The dimensional |]rT| and 
Pauli-Villars ^| regularizations are typical examples of this type. 

We begin with the dimensional regularization. The idea is that the volume of (|5.5| ) 
depends on the number of dimensions D — 1 and it diverges as e^~^ where e is the cutoff 
at small p {p > e). Formally, ii D < 2, the integral at small p converges so one can use D 
as a regularization parameter. Consider as an example a massive scalar fielcQ described 
by Eq. ( p.l8|) with V = m?. Near the horizon the diagonal element of the trace of the 
heat kernel of (B.4|) behaves as ITT 



-HH 



diag (47rt)(^-i)/2 



1 + M + M^ + --) 



1 



6i~«:V(^--e)i?, h = 0{p'). 



(5.6) 



(5.7) 



Other terms in hi do not result in singularities at D ^ 4. The integration in the trace 
near p = can be written as 



where V is the curvature invariant computed on E 



l + ^p^P 



(5. 



V = 2P^"'P^PR^xup - P^'^R 



<^fiu- 



(5.9) 



In case of stationary but not static space-times V has to be computed by using the 
curvatures of the corresponding fiducial static background Ai. However, as we have 



mentioned, curvatures P on 7W and 7V1 coincide ^8 
Equations (K 



]q) can be used to derive the divergent part of the trace 



Tr e 



-hH 



r(i 



m 



D-4 



div 



(47r)(^-i)/2 2fi:t3/2 



22 



The mass term elim.inates divergencies at large p. 



m 



l-^ 



R]t 



V 



(5.10) 



38 



From ( |5.10| ) one can also get the spectral density of H^ in dimensional regularization (see 

(EUD) 



'^div(^'^) 



r(i-f) ^D-A 



(47r)^/2 
and from it the free energy 



2 im' 



^]R 






V 



(5.11) 



^div(/^'^) 



r 1-4^ 



2) TT-m 



2^D-4 



(47r)^/2 3k(3^ 



m 






V 



-,-^)R 



(5,12) 



The divergencies in (|5.11|) , (|5.12|) appear as poles at D = 4. This is a typical property of 
dimensional regularization. 

Expression (|5.11|) is suitable to go to the Pauli-Villars (PV) regularization. Let us 
introduce 5 additional auxiliary scalar fields (P V partners) : 2 with masses M^ and 3 with 
masses M^, and impose two restrictions 



where 



/(l) = /(2) = 0, 



/(p) = m2^ + E^f-E(^;r^ = o 



(5.13) 



(5.14) 



Equation ( |5.13| ) has the solutions Mi^2 = A/3/?~+m^, M[ 2 = V/x^ + m^, M3 = A/4/?~+7n^ 
where /i is a regularization parameter, see [^. The regularized density of states in PV 
method is defined as 



$(0;, /i) = $(^, m) + E -^l^, Mk) - E "^(^^ ^r) 



(5.15) 



where $(u;,Mfc) and $(a;,M^) are spectral densities of PV partners. Fields with masses 
Ml give negative contributions to the regularized density (as if these scalars had a wrong, 
Fermi, statistics). 

Suppose now that we start with dimensionally regularized theory. At D 7^ 4 the 
spectral density for each PV partner is a finite quantity given by (|5.11| ). Although each 
of these densities has a pole the poles are cancelled in the total density ( |5.15| ) by virtue 
of the constraints ( |5.13| ). Hence, one can put in ( |5.15| ) D = A and get in the limit fi ^ m 



%v(^'/^) 



1 



(47r)2K Js 






(5.16) 



where c = ln(729/256). The corresponding expression for the divergent part of canonical 
free energy is 



^div(/^'/^) 



48fi;/?2 7s 



/i- 



2c/i^ + In ^ 



1 47r2 



im? V 60 K^/S 



V+ [--i]R~m 



(5.17) 
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The Pauli-Villars regularization effectively results in a cutoff of the integrals near the 
horizon at the proper distance comparable to the inverse mass of PV partners. This fact 
allows the following interpretation. Near the horizon, where the local temperature be- 
comes greater than /i, PV fields become thermally excited and by virtue of the constraints 
( |5.13|) their contribution exactly cancels the contribution of the physical field. 

Usually dimensional and Pauli-Villars regularizations are used to regularize the ul- 
traviolet divergencies. As we have seen, in case of black holes they can be also used to 
avoid infrared singularities near the horizon. We will show now that this is not an acci- 
dent and the divergences of canonical free energy in these regularizations match precisely 
ultraviolet divergences of Euclidean free energy. 

5.4 Heat kernel expansion and conical singularities 

Let us denote by Aia a Euclidean manifold which possesses a one-parameter group of 
isometries. Suppose that S is a fixed-point set of these isometries and that J^^ has a 
structure (7^/3 x i?-^"^ near S. Here C^ is a two-dimensional conical space (the value 
a = 2-11 corresponds to a plane). 

....1 .,..„.„... ^f „ T „.!„„ .„... .. A.,^ 



One has the following heat kernel expansion of a Laplace operator on A^f 



Tr e-*^- ~ j—^^ i:(A„ + A^,„)t". (5.18) 

Here An are standard coefficients defined on the regular domain of A4^. Additions A^ „ 
appear due to conical singularities and accumulate information about local geometry near 
S. Ai3^n are local functionals on S expressed in terms of powers of the Riemann tensor 
and its derivatives projected with the help of P^'^. The quantities PAp^n are polynomial 
in /3~^. For the operator Le = — V^ + ^R one has 

A^,o = l , A^,i = T^(7'-l)i, (5.19) 
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^^■■'=^.LU''-'^^^^''-'^{6-¥. 



(5.20) 



TT 
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where 7 = 27r/(K/5), A is the surface area of Euclidean horizon S and V is the invariant on 
S defined by ( |5^) . Laplacians on Ca and the coefficient Ajs^i were discussed by Cheeger 
|63[| , (see also |^,||3l|). The heat kernel expansion on spaces with fixed-point sets was 
analyzed by Donnely ||6^. Heat kernel expansion (|5.18|) on Ai^ has been studied in Refs. 



3^ , ||65[| , ||66[| where the explicit form of ( ^.20|) was found. A summary of results for 



asymptotics for higher spin fields is presented in Ref. ||5^ . 
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5.5 Infrared /ultraviolet duality 

Now we define the one-loop effective action W^{l3) on A^f by using dimensional regular 



ization, I.e. 



1 /-OO rif ,, 



W^{(3) = -\ndet{LE + m^) = - —e"*'" Tr e"*^^. (5.21) 

2 Jot 



For D y^ 4 the singular part of (|5.21| ) can be found from (|5.18 



1 /-OO rlf p-tm-' 

'*'div('3) = -5i JJi^^ABo + tB^^fB,), (5.22) 

where Bk = Ak + Ap^k- On a regular Euclidean manifold the ultraviolet divergences are 
determined only by Aq, Ai, ^2 which are proportional to the period (3. So they appear 
only as divergencies E^ ^^j^ of the vacuum part Eq of F^{P), see ( |1.16| ). In case of conical 
singularities the ultraviolet divergences are polynomials in P~^ because of A^ ^ and they 
also appear in the thermal part of F^{(3). The divergent part F^- {(3,D) can be found 
from ( ^:i9D , (|]2g) and (|]22|). Note that F|. (/?,£)) is a local functional determined 



by curvatures and geometrical invariants {A, R, V) at the Euclidean horizon. For this 
local functional the Wick rotation described in section 3.1, Eqs. (|3.2| ), has a well-defined 
meaning and it can be used to go in F^- (/5, D) from Euclidean to Lorenztian theoryQ. As 
one can further show, after the Wick rotation the thermal part of F^. (/5, D) is identical 
to the divergent part of the canonical energy ( |5.12| ). This coincidence also takes place in 
the Pauli-Villars regularization, so one can write 



^div(/^' ^) - ^o,div('^) = ^div(/^' '^)' (5-23) 

where 5 is a regularization parameter [6 = D — 4 for dimensional regularization and 
6 = /i~^ for PV regularization). 

Relation ( p.23| ) is remarkable. It demonstrates a duality between large-volume in- 



frared divergences which appear in F^{(3) and ultraviolet-type divergences which appear 
in F^{(3) as a result of curvature singularities on the Euclidean horizon. This duality is 
manifest in the Pauli-Villars regularization where the inverse scale of Pauli-Villars masses 
/i~^ plays the same role as the cutoff parameter of integrals near the horizon. 

We conclude our discussion with very brief comments on the problem of statistical 
explanation of the Bekenstein-Hawking entropy S^^ of black holes. Note that the entropy 
of a field near the black hole horizon can be computed by the standard formula S = 
(3'^dpF'~^ at the Hawking temperature (3~^ = k/{2t[). By using ( p. 17 ) one can show 



that to leading order S ~ A/ 5"^ where A is the surface area of the black hole horizon 
and 5 is the cutoff. The infrared/ultraviolet duality enables one to identify 5~^ with 
a high-energy cutoff and, in particular, with the Planck mass. This fact is interesting 

^■^This procedure is analyzed in some detail in |^ for a Kerr-Newman black hole. 
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because it can be used as one of possible explanations of the Bekenstein-Hawking entropy 
S^^ = A/{4:G). Indeed, suppose that Einstein gravity is entirely induced by quantum 
effects of some fundamental underlying theory. For simplicity one can assume that the 
degrees of freedom (constituents) of this theory are some field variables. Then the induced 
Newton constant is G ~ 5^ where 6~^ and S^^ ~ A/S"^. Because of the duahty property 
the same cutoff 6 which determines the induced Newton constant appears also in the 
entropy of constituents S and automatically S has the same order of magnitude as S^^ . 



A similar mechanism may be realized in open string theory |68]. Further discussion of the 
black hole entropy problem in induced gravity and concrete realizations of this idea can 
be found in Ref. [^ . 



6 Resume 

The aim of these notes was to overview a number of new results and methods in finite- 
temperature field theory in classical external stationary backgrounds. We have also dis- 
cussed their physical applications, including quantum effects near rotating black holes. 
The physics here is reduced to a class of spectral problems which depend on the spec- 
tral parameter in a non-linear way. This also represents an interest form mathematical 
point of view because some results of the spectral theory such as heat kernel asymptotics 
can be extended to this class of problems. One of our motivations was to find a con- 
nection between canonical and Euclidean methods in statistical mechanics in situations 
more non-trivial than static external fields. It is interesting to note that the Euclidean 
gravity yields an alternative derivation of the spectral asymptotics. This enables one to 
avoid a complicated analysis of inner products and measures which has to be done in the 
Lorenztian theory. This fact hints that our results might be extended to a larger class of 
spectral problems. 
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